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We develop a model for describing excited mesons decay into three mesons. The properties of the 
excited mesons can be extracted with this model. The model maintains the three-body unitarity 
that has been missed in previous data analyses based on the conventional isobar models. We 
study an importance of the three-body unitarity in extracting hadron properties from data. For 
this purpose, we use the unitary and isobar models to analyze the same pseudo data of jp — > 
n + n + n~n, and extract the properties of excited mesons. We find a significant difference between 
the unitary and isobar models in the extracted properties of excited mesons, such as the mass, 
width and coupling strength to decay channels. 
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1. Introduction 

Hadron properties such as quantum numbers (spin, parity, etc.), mass and (partial) width have 
been long studied as a subject called hadron spectroscopy. The hadron properties provide important 
information for understanding internal structure of the hadron and dynamics which governs it. The 
dynamics here is of course QCD in its nonperturbative regime. The hadron properties can be 
extracted from data through a careful analysis, in many cases, partial wave analysis (PWA). Thus 
it is essential for hadron spectroscopy to have a reliable theoretical analysis tool. 

Here, we are interested in analyzing data in which excited mesons decay into three mesons. 
Such an analysis has been done, for example, in Ref. [jl|] in which kN — > M*N — > TinnN (M*\ an 
intermediate excited meson) is analyzed to extract the properties of M*. The analysis tool which 
has been conventionally used for this kind of analysis is the so-called isobar model. In the isobar 
model, it is assumed that two mesons form a resonance (fo,p,fi, etc.), and the third meson is 
treated as a spectator, as illustrated in Fig. |j] (a). A naive question that arises here is how significant 
the missing three-body unitarity and coupled-channels effect are. In order to take care of the three- 
body unitarity, we need to consider the Z-diagrams [Fig. [j] (b)]. Resumming all-order multiple 
scattering diagrams due to the Z-diagrams, as shown in Fig. 0, we obtain a M* decay amplitude 
that maintains the three-body unitarity, and also coupled-channels effects as a consequence. 

In this contribution, we would like to address a question how the three-body unitarity makes 
a difference in extracting hadron properties from data [||. To address this question, we construct 
both unitary and isobar models, and fit them to the same data. Then we extract and compare M* 
properties from the two models. We will conduct this analysis for the yp — > 7i + 7i + n~n reaction. 
The CLAS Collaboration tried to find an exotic meson in this reaction using an isobar model [Q], 
and the GlueX Collaboration plans to do a more elaborate study on this [Q]. To pin down the 
existence of exotic states, a reliable analysis tool is essential. Therefore, it is interesting to study 
differences between the unitary and isobar models for this reaction. 

The organization of the rest of this report is as follows: We discuss our unitary model in Sec. ^. 
We also define an isobar model in the same section. In Sec. f|, we determine parameters contained 
in the unitary model using a result from a 3 Po model calculation [||]. Then we use the unitary model 
to generate pseudo data (Dalitz plots) for yp — > 7t + 7i + 7i~n, and fit the data with the isobar model. 
M* properties are extracted from the two models, and are compared. 

2. Unitary coupled-channels model 

We assume that the yp — > 7t + 7i + 7i~n reaction proceeds through yp^M*n andM* — > 7i + 7i + n~, 
schematically shown in Fig. || (left). There are some other mechanisms, such as an excitation of 
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Figure 1: M*-decay amplitude. 
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Figure 2: (Left) Schematic yp — > K + K + K~n mechanism considered in this work; (Middle, Right) M* decay 
amplitude calculated with the unitary model. The label "t" is the nR amplitude, and the gray circles near R 
stand for intermediate loops of KK and KK. The isobar model considers only the middle diagram. 

nucleon resonances and/or the Deck effect, that lead to the same final state. We do not consider 
these backgrounds because we do not deal with actual experimental data as we mentioned. For our 
limited purpose of investigating effects of three-body unitarity, we expect that these backgrounds 
would not play an important role. Thus in the following, we give a discussion on the M* — > 3k 
amplitude described with our unitary model (a detailed discussion is given in Ref. |0]), and then 
discuss M* production mechanism and its propagation. The yp — >■ K + K + K~n amplitudes are prod- 
uct of these ingredients, and are used to calculate Dalitz plots distributions of the three pions. 

M* decay amplitude We assume that an M* decays through M* — > kR — » kkk (Fig. [T]) where R 
stands for a bare state formed by two pions. Two pions from /?-decay can form R again. Through 
this type of rescatterings, R is dressed by the pion cloud to develops a resonance pole corresponding 
to, for example, fo (600), p (770), etc. Therefore, we include these resonances in our model as poles 
in the KK scattering amplitude; we do not use Breit-Wigner functions that has been often used in 
previous isobar-model analyses. Meanwhile, one of two pions from /?-decay can also interact with 
the other pion to form R. We call this mechanism the Z-diagram [Fig. |T](b)]. Thus a basic ingredient 
in our model for describing the final state interaction of M* -decay is the R^ KK interaction. This 
interaction can be fixed by analyzing KK scattering data. We developed a simple coupled-channels 
(KK,KK) model for the KK scattering For each partial wave, KK (or KK) — > KK (or KK) 
potentials are given by ^-channel exchange of bare R states, and the partial wave amplitude is 
obtained by solving three-dimensional Lippmann-Schwinger equation. We obtained a reasonable 
description of empirical KK phase shifts and inelasticities for (L,I) = (0,0), (1,1), (2,0) partial 
waves (L: orbital angular momentum; /: total isospin) from the threshold to 2 GeV. We extracted 
pole positions from the partial wave amplitudes, and found that they are in good agreement with 
the listing in the Particle Data Group. 

Having fixed the R -H- KK vertices, we now have a coupled-channels scattering equation for 
kR amplitudes, and the Z-diagrams work as the driving force. Different kR channels are coupled 
via the Z-diagrams. The R Green function contains self energies from the KK and KK loops. Both 
Z-diagrams and kR Green functions contain three-pion unitarity cut, and the resultant three-pion 
amplitude (from the kR amplitude) satisfies the three-body unitarity. The three-pion unitarity cut in 
the Z-diagrams makes it difficult to solve the kR scattering equation with the standard subtraction 
technique. We handle the problem with the Spline method which is discussed in detail in Ref. [^J. 

The kR amplitudes are used to calculate the M* decay amplitudes, as seen in Fig. |2|(middle,right), 
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by the convolution with M* — > kR vertices. By retaining only the middle diagram of Fig. ^| the M* 
decay amplitude is similar to those used in the isobar models. Thus we define our isobar model by 
considering only this mechanism. 

Green function and production mechanism of M* The M* Green function for the unitary 
model is given by G~\W) = W - m° M , - L(W), where W is the total energy of the three pions 
in their CM frame, m^* is the bare mass of M*, and £(VK) is the self energy. The self energy 
contains the kR amplitudes, so that the self energy has a consistency with the M* decay amplitude 
(see Ref. [Q] for expressions). The three-body unitarity requires this consistency. Pole positions 
corresponding to M* resonances are a quantity of interest. We search for pole positions (Mr) 
that satisfies G~ 1 (Mr) = 0. The pole search needs an analytic continuation of the amplitudes to 
the complex energy plane. A method of the analytic continuation suitable for our unitary model 
has been developed in Ref. [gj. For the isobar model, on the other hand, we use a Breit-Wigner 
function, as has been done in most isobar-model analyses. 

For the M* production mechanism, we assume a ^-channel pion exchange, and the pion inter- 
acts with a p-meson from the incident photon, invoking the vector-meson dominance. Although 
this is a rather simple assumption and often not very realistic, this model still works enough for our 
purpose of studying effects of the three-body unitarity. We use this M* production mechanism for 
both of the unitary and isobar models. 

3. Numerical results 

Having described our unitary and isobar models, we now use them for a numerical analysis. 
As mentioned in the introduction, we want to learn how three-body unitarity makes a difference 
in extracting M* properties from Dalitz plot analysis. For this purpose, we take the following 
procedure: (i) We determine M* — > kR couplings for the unitaiy model using a prediction from the 
3 Po model for partial width for M* — > nR decays [^] ; (ii) We generate pseudo data with the unitary 
model; (iii) We fit the data with the isobar model; (iv) We extract the M* properties from the two 
models. Details are discussed in the following paragraphs. 

We consider partial waves and M*s that have been included in the CLAS analysis of yp — > 
K + n+n-n%. They are J PC = 1++ [a { (1230), a x (1700)], 2++ [a 2 (1320), a 2 (1700)], 2"+ [7T 2 (1670), 
7T 2 (1800)], and 1~+ [^(1600)]. The 3 P model, a quark model, can predict partial widths for 
M* — > 7iR decays We use the prediction to fix M* —> nR coupling constants, assuming that 
the phases are +1. With the parameters fixed in this way, we run the unitary model to generate 
Dalitz plots. We chose the following kinematics: E y =5 GeV (Ey\ incident photon energy in the 
laboratory frame); t = —0.4 GeV 2 (?: squared four-momentum transfer of the nucleon); 0.8 GeV 
< W < 2 GeV. Finally, three-pion orientation is chosen so that the plane formed by the three pions 
is perpendicular to the final nucleon momentum in the three-pion CM frame. In Fig. g, we show 
a spectrum obtained by integrating the three pion distribution within the specified kinematics. As 
seen in the figure, different M*s play major roles at different W regions, depending on their masses 
and widths. At each W, the final three pions are distributed in a pattern characteristic to the as- 
sociated M*. For example, in Fig. ^| (left), we show a Dalitz plot at W =1 GeV where the main 
contribution is from the ai(1230) decay. A three-dimensional view of the same Dalitz plot is also 
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Table 1: M* masses and widths from the unitary (UT) and the isobar (IB) models. 





ai(1230) 


a 2 (1320) 


7Tl(1600) 


tt 2 (1670) 


ai(1700) 




UT IB 


UT IB 


UT IB 


UT IB 


UT IB 


Mass (MeV) 


937 1096 


1263 1269 


1599 1599 


1784 1798 


1658 1687 


Width (MeV) 


496 684 


42 56 


8 8 


456 502 


106 99 



shown by the solid curve in Fig. |] (right). On the same figure, the dotted curve is obtained by 
simply turning off the Z-diagrams from the M* decay amplitudes. Thus the difference between the 
solid and dotted curves is the effect of the Z-diagrams (of the M* decay amplitude) on the Dalitz 
plot. The Dalitz plot distributions are fitted with the isobar model. Fitting parameters are (complex) 
M* TiR couplings and cutoffs, the Breit-Wigner masses and widths in the M* Green functions. 
The Breit-Wigner width also contains an adjustable parameter to fit the energy dependence. We 
also include a flat (in three-pion distribution), non-interfering background term, as usually done 
in the isobar-model analyses. The strength of the flat background is also a fitting parameter. The 
quality of the fit can be seen by comparing the dashed curve (fit) with the solid curve on the right 
panel of Fig. ^. The isobar model reasonably fits the Dalitz plots from the unitary model. 

We now discuss the M* properties extracted from the unitary and isobar models. The mass 
(width) of anM* of the unitary model is defined with the pole position, Mr, by Re[M#] (—2 Im[M#]). 
These are compared with Breit-Wigner masses and widths from the isobar model in Table [I]. While 
we can see a good agreement for narrow M* [a2(1320), n\ (1600)], there are also significant differ- 
ences for broader M* [ai(1230), 7T2(1670)]. Another interesting quantity is the coupling strength 
of M* to each decay channel. Our result for a\ (1230) — > TiR vertices are shown in Fig. |] as a func- 
tion of the pion momentum, and the unitary and isobar models are compared. The difference is 
significant in both the strength and the phase. These results indicate the importance of considering 
the three-body unitarity in extracting the M* properties. 




Figure 3: W-dependence of integrated Dalitz plots from the unitary model. The M* — > nR couplings in the 
model have been fixed using the 3 Pq model prediction. M's contributing to the peaks are indicated. 
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Figure 4: (Left) Dalitz plot distribution of 7T+7T+7Z; in the yp — > n + n + n n reaction at W = 1 GeV. (Right) 
3 -dimensional view of the Dalitz plot distribution shown on the left. 
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Figure 5: The strength (phase) of M* — > nR vertices on the upper (lower) half. The unitary (isobar) model 
is displayed by the solid (dashed) lines. (Left) a\ (1230) — > Kfy vertex in / = 1 (/: orbital angular momentum 
between % and R). (Right) a\ (1230) — > Tip vertex in I = 0. 
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